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Description 

A GEOPHYSICAL METHOD FOR 
QUANTIFICATION OF DENSE 
NON-AQUEOUS PHASE LIQUIDS IN THE 

SUBSURFACE 

Cross Reference to Related Applications 

[0001] This disclosure claims the benefit of pending provisional 
patent application Serial No. 60/319,380, filed on July 5, 
2002. 

Background of Invention 
[0002] FIELD OF INVENTION 

[0003] This invention relates to the quantification of contami- 
nants in soil and more particularly to a geophysical 
method for quantifying the volume of dense non-aqueous 
phase liquids in the subsurface. 

[0004] BACKGROUND OF THE INVENTION 

[0005] Dense non-aqueous phase liquids (DNAPL) are common 



and highly problematic contaminants of concern in 
ground water. They are of concern at dissolved concentra- 
tions orders of magnitude lower than their solubility, they 
can migrate many meters vertically in soils within hours or 
days, they pool at horizontal porosity and permeability 
boundaries, and their migration is little effected by 
groundwater flow or gradients. 
[0006] Quantifying the volume of non-aqueous phase liquids in 
the subsurface is difficult. To date, methods have been 
developed which can detect DNAPL in the subsurface but 
these methods are not effective in quantifying the amount 
of DNAPL present. As cleanup costs at DNAPL sites typi- 
cally run in the millions of dollars, accurately quantifying 
the volume of DNAPL present is very important. There is a 
need for a method to locate and quantify the extent of or- 
ganic contaminants in the non-aqueous phase in order to 
develop strategies to remediate the DNAPL sources. Site 
remediation is dependent upon reducing a certain per- 
centage of the DNAPL contaminant mass that is the source 
for the dissolved contaminant plume. In order to estimate 
the amount of mass that needs to be removed, the per- 
cent saturation of DNAPL within the subsurface needs to 
be reasonably quantified. 



[0007] Accordingly, what is needed in tlie art is a metliod tliat 
can directly detect and quantify the volume of DNAPL 
present in a saturated soil. 

[0008] It therefore, to the effective resolution of the afore- 
mentioned problems and shortcomings of the prior art 
that the present invention is directed. 

[0009] However, in view of the prior art in at the time the present 
invention was made, it was not obvious to those of ordi- 
nary skill in the pertinent art how the identified needs 

could be fulfilled. 
Summary of Invention 

[0010] jhe present invention is a method of detecting and quan- 
tifying non-aqueous phase liquids in a soil sample. In a 
preferred embodiment, the method includes the steps of 
measuring the bulk dielectric constant of the sample, 
measuring the bulk density of the sample, estimating the 
bulk porosity of the sample from the measured bulk den- 
sity, and calculating the non-aqueous phase liquid satu- 
ration from the measured bulk dielectric constant and the 
estimated bulk porosity. 

[0011] In an additional embodiment, 100% water saturation and a 
matrix density 2.65g/cc are assumed when estimating the 
bulk porosity from the bulk density measurements of the 



sample. 

[0012] In another embodiment, the non-aqueous phase liquid is 
a dense non-aqueous phase liquid. 

[0013] In yet another embodiment, the soil sample is in the sub- 
surface. The steps of measuring the bulk dielectric con- 
stant and measuring the bulk density are performed in 
situ and the bulk dielectric constant of the subsurface is 
measured at a plurality of depths and bulk porosity of the 
subsurface is measured at the same plurality of depths 
and the bulk dielectric constant measurements at the plu- 
rality of depths and the bulk porosity measurements at 
the same plurality of depths are used to calculate the 
non-aqueous phase liquid saturation at the plurality of 
depths. 

[0014] In a preferred embodiment of the present invention, the 
sample is in the subsurface and the measurements are 
performed in situ, whereby a logging well and an active 
gamma logging tool are provided to measure the bulk 
density of the subsurface at a plurality of vertical intervals 
along the length of the logging well. A cone-penetrometer 
tool (CPT) is provided having dielectric constant measure- 
ment capability. The CPT is used to measure the dielectric 
constant of the subsurface at a second plurality of vertical 



intervals. The bulk density measurements and the dielec- 
tric constant measurements are interpolated to establish a 
series of measurements at a plurality of consistent vertical 
intervals. The bulk porosity of the subsurface at the plu- 
rality of consistent vertical intervals is then estimated 
from the interpolated bulk density, assuming 100% water 
saturation and a matrix density of 2.65g/cc. The percent 
non-aqueous phase liquid saturation is then determined 
at the plurality of consistent vertical intervals from the in- 
terpolated dielectric constant and the interpolated bulk 
porosity. 

[0015] In yet another embodiment, a cone-penetrometer tool is 
provided. The cone-penetrometer tool includes a bulk 
density measurement device and a dielectric constant 
measurement device. With the bulk density measurement 
device and the dielectric constant measurement device 
existing on the cone-penetrometer tool the bulk density 
measurements and the dielectric constant measurements 
of the subsurface can be obtained simultaneously at a 
plurality of vertical intervals within the subsurface. Inter- 
polation is not necessary and the measurements can be 
used to estimate the bulk porosity of the subsurface and 
for the subsequent calculation of the percent non- 



aqueous liquid saturation of the subsurface at the plural- 
ity of vertical intervals. 

[0016] It is to be understood that both the foregoing general de- 
scription and the following detailed description are ex- 
planatory and are not restrictive of the invention as 
claimed. The accompanying drawings, which are incorpo- 
rated in and constitute part of the specification, illustrate 
embodiments of the present invention and together with 
the general description, serve to explain principles of the 
present invention. 

[0017] These and other important objects, advantages, and fea- 
tures of the invention will become clear as this description 
proceeds. 

[0018] The invention accordingly comprises the features of con- 
struction, combination of elements, and arrangement of 
parts that will be exemplified in the description set forth 
hereinafter and the scope of the invention will be indi- 
cated in the claims. 
Brief Description of Drawings 

[0019] For a fuller understanding of the nature and objects of the 
invention, reference should be made to the following de- 
tailed description, taken in connection with the accompa- 
nying drawings, in which: 



[0020] FIG. 1 is a table of dielectric constants for selected mate- 
rials. 

[0021] FIG. 2 is an equation used to calculate the porosity as- 
suming 100% water saturation from the results of the bulk 
density measurements. 

[0022] FIG. 3 is an equation used to calculate formation bulk 
density from porosity, the matrix density of quartz and 
the average density of pore fluid. 

[0023] FIG. 4 illustrates the results of the bulk density calcula- 
tions using the equation of Fig. 3. 

[0024] FIG. 5 is an equation used to calculate the percentage of 
change in bulk density over a 100% change in DNAPL sat- 
uration. 

[0025] FIG. 6 is an equation used to calculate bulk dielectric. 

[0026] FIG. 7 is an equation showing the relationship between 
chlorine concentration and the court rate for chlorine. 
[0027] FIG. 8 shows dielectric response to porosity. 

[0028] FIG. 9 shows dielectric response to DNAPL. 

[0029] FIG. 10 illustrates the results of a natural gamma log used 

to determine the site stratigraphy. 
[0030] FIG. 11 is a graph illustrating the bulk dielectric constant 

relative to elevation of the site. 



[0031] FIG. 12 is a graph illustrating the porosity relative to ele- 
vation of the site. 
[0032] FIG. 13 is a graph illustrating the neutron log for the well. 

[0033] FIG. 14 illustrates the calculated dielectric constant rela- 
tive to the field values of porosity. 

[0034] FIG. 15 illustrates the difference between the calculated 
and observed values of the dielectric constant. 

[0035] FIG. 16 illustrates the quantitative measure of the volume 
of DNAPL present as a fraction of the total volume. 

[0036] FIG. 17 is an illustrative plot of DNAPL saturation versus 
depth. 

[0037] FIG. 18 is a flow diagram illustration the method for de- 
tection and quantification of a non-aqueous phase liquid 

as described by the present invention. 
Detailed Description 

[0038] The present invention utilizes the measurements of the 
bulk dielectric and the bulk density of the subsurface to 
detect and quantify the volume of DNAPL present at a 
sample point. 

[0039] The dielectric constant of a medium is the ratio squared 
of the electromagnetic propagation velocity in a vacuum 
relative to the velocity in that medium. The variables 



which can affect this electrical property include soil tex- 
ture, water content and density. The dielectric constant 
values that can be used to identify materials in the sub- 
surface are presented in Table 1 of Fig. 1. It is known that 
DNAPLs have very low dielectric constant values (3-4) 
compared with water (80). Because the dielectric values 
vary greatly between water and DNAPL, displacement of 
pore water with DNAPL significantly changes the mea- 
sured value of the bulk dielectric constant. Replacing wa- 
ter in pore spaces with DNAPL will cause a decrease in the 
dielectric constant. 
[0040] An electronic piezocone (cone penetrometer, or CPT) is 

used to provide hydrogeologic profiling for environmental 
and geotechnical projects. The CPT provides a rapid, reli- 
able and economical means of determining soil stratigra- 
phy, relative density, strength and hydrogeologic informa- 
tion in contaminated soils, without generating contami- 
nated waste cuttings. The present invention utilizes an 
electronic piezocone or cone penetration test (CPT) to 
measure variations in dielectric constant with depth. The 
CPT is employed along with stratigraphy and bulk resistiv- 
ity. The numbers of CPT holes are chosen to cover as 
much of the area extent of the suspected source area as 



possible and to penetrate at least to the top of the first 
significant low permeability layer. The spatial arrangement 
of the CPT holes within the area is based on site layout 
and accessibility for both the CPT equipment and drilling 
rigs. Gamma-gamma logging, or bulk density, logging is 
employed to measure the bulk density of the sample. An 
active gamma source generates high energy gamma radi- 
ation, some of which is reflected back to the logging tool 
by the surrounding soil. The amount of back-scattered 
gamma radiation is proportional to the bulk density of the 
soil. The gamma-gamma logging is performed in moni- 
toring wells, which are preferred to be next to a CPT hole. 
Current technology requires density logging in boreholes 
adjacent to the CPT locations, but the dielectric and den- 
sity logging tools could be placed on the CPT tool, elimi- 
nating the need for the boreholes. 
[0041] It is known that the dielectric constant varies with both 

porosity and DNAPL saturation. A decrease in porosity will 
also cause a decrease in the dielectric constant, so to use 
the dielectric constant to quantify the volume of DNAPL 
present; it is necessary to also measure porosity. The 
present invention addresses this problem by using the 
gamma-gamma log. In order for this solution to be feasi- 



ble, the CPT holes need to be located next to monitoring 
wells that will be used for the gamma-gamma logging. 
The gamma-gamma log uses a sonde that is lowered 
within the borehole until it is in contact with the formation 
wall or well bore. A gamma-ray source attached to the 
probe emits gamma rays that pass into the formation and 
bombard the electrons within the formation. Gamma rays 
that are scattered back to the detectors are counted and 
used to calculate the bulk density. Using the known den- 
sities of water and quartz, porosity can be calculated from 
bulk density. Once variations in porosity are identified, 
variations in dielectric constant can be correlated to 
changes in DNAPL. The results of the bulk density mea- 
surements are used to calculate the porosity assuming 
100% water saturation. This equation is provided in Fig. 2. 
[0042] A complication encountered is that the presence of DNAPL 
changes the bulk density, and assuming 100% water satu- 
ration leads to an error in the calculation of porosity. 
However, the effect of DNAPL on bulk density is much less 
than on the bulk dielectric. The percent difference in bulk 
density values from variations in DNAPL saturation is very 
small and can be illustrated using the equation of Fig. 3, 
assuming 100% pore saturation. As an example, using the 



equation of Fig. 3, an estimated porosity of 0.2, and per- 
centages of DNAPL of 0%, 50% and 100%, the results for 
changes in bull< density are provided in Fig. 4. The per- 
cent change in bull< density (maximum change) over a 
100% change in DNAPL saturation can be calculated as 
shown in Fig. 5. It has been shown that DNAPL can occupy 
40-70% of the bulk pores space. If the maximum possible 
DNAPL saturation is approximately 70%, the percent error 
in calculating the bulk densities as a result of the pres- 
ence of DNAPL will range from 0-2.8%. 
[0043] In order to determine how the bulk dielectric will change 
with variations in DNAPL saturation, a mixing model 
equation is applied. This equation calculates bulk dielec- 
tric constant in a three-phase system assuming 100% 
pore saturation and is dependent upon, porosity, medium 
geometry, and the dielectric constants of water, soil and 
DNAPL. The equation used to calculate the bulk dielectric 
constant for a DNAPL, water, and soil system is obtained 
by substituting a term for DNAPL for the term for air in 
the mixing equation for a water, soil, and air system as 
provided in Fig. 6. Using the equation of Fig. 6, the bulk 
dielectric can be calculated using a range of porosity val- 
ues for different percent DNAPL saturation values. For a 



porosity of 0.20, the bulk dielectric constant decreases by 
about 35% as the DNAPL saturation (% of available pore 
space) increases from 10% to 70%. These calculations il- 
lustrate that DNAPL saturation has a significant effect on 
the bulk dielectric. The bulk dielectric changes by 80% 
from varying percent DNAPL saturation from 0 to 100% 
DNAPL at a constant porosity. The percent change in bulk 
dielectric with varying porosity at three selected DNAPL 
saturation values ranges from 0 (porosity = 0) to a maxi- 
mum of 219 (porosity = 0.7). There is a direct relationship 
between bulk dielectric and porosity while there is an in- 
verse relationship between bulk dielectric and percent 
DNAPL saturation. 
[0044] The previous calculations with respect to the equation of 
Fig. 3, demonstrate the small effect that DNAPL saturation 
has on bulk density values. The range in percent differ- 
ences of bulk density values were from 0 - 2.8% over a 
range of percent DNAPL saturation values from 0-70%. 
The presence of DNAPL will not change the bulk density 
significantly enough to have a noticeable effect on the es- 
timated porosity. Therefore the estimated porosity ob- 
tained from the gamma-gamma log results can be used in 
Fig. 6 with the bulk dielectric values to determine the 



DNAPL saturation. 

[0045] The fractional DNAPL volume was calculated using the 
equation of Fig. 6, the observed bulk dielectric constant, 
porosity, and the dielectric constants for water, quartz 
sand, and DNAPL The dissolved phase constituents of the 
sample are used to determine the dielectric constant for 
DNAPL. Porosity, water content and DNAPL content are 
expressed as fractional volumes, which are the volume of 
the pores, the water or the DNAPL divided by the total 
volume. DNAPL saturation refers to the percentage of 
available pore space occupied by DNAPL. For example, a 
fractional DNAPL volume of 0.20 in a soil with a porosity 
of 0.40 yields a DNAPL saturation of 50%. 

[0046] The present invention provides a non-intrusive, quantifi- 
able and defensible calculation of the percent of DNAPL 
saturation within the subsurface. Applied spatially, the re- 
sults can be used to estimate the amount of DNAPL mass 
present in the subsurface and could be used in a mass 
transport model to predict how the dissolved plume will 
change over time based upon a certain percentage of 
mass removal. Knowledge of the amount of DNAPL 
present will also aid in choosing the appropriate remedia- 
tion technology. The practice can be applied at DNAPL 



contaminated sites and a gamma-gamma logging tool 
could be installed with the CPT and resistivity probe to 
combine the technologies, thus eliminating the depen- 
dence on monitoring wells and the costs associated with 
their installation. 

[0047] The following example is provided as exemplary of the 
method of the present invention and is not intended to 
limit the scope of the inventions or claims thereto. 

[0048] The method of the present invention was employed for 

the quantification of DNAPL in saturated soils at a test site 
known to be heavily contaminated with DNAPL, principally 
TCE and DCE. The site was located on a large barrier is- 
land complex in the Southeastern U.S. The presence of 
DNAPL was confirmed by water-quality sampling using di- 
rect push tools to define the study area as shown in Fig. 7 
where the contaminated areas 1 and clear locations 5 are 
identified. High concentrations, greater than 100,000 \xg/ 
I, of dissolved phase TCE occur at the site as shown in the 
table of Fig. 8. Dissolved concentrations in excess of 1% of 
saturation are considered to indicate the presence of 
residual DNAPL. 

[0049] The general site stratigraphy consists of marine sands and 
silts with some clay. Land surface elevations are about 5 



to 8 ft above sea level (msl). Clean, well-sorted sands are 
present from the surface to elevations of 0 to +2 ft msl. 
Three or four silty layers a few feet thick each are present 
in the upper 35 to 40 ft, and a persistent clayey, silty sand 
at elevations of 40 to 45 ft msl forms the bottom of the 
upper surficial aquifer. The water table is at an elevation 
of about +2 ft msl under the site. 

[0050] The original source of the DNAPL, and the amounts dis- 
posed of, are not known, but the dissolved concentrations 
suggest that the residual DNAPL is confined to a roughly 
oval area about 100 ft wide and 200 ft long. Highest dis- 
solved TCE concentrations are found between elevations 
of 7 to 22 ft msl, with concentrations generally falling be- 
low 100,000 ng/l at elevations below 30 to 35 ft msl. 

[0051] Five wells were constructed to obtain the borehole logs. 
The wells were drilled using the rotosonic method to min- 
imize soil disturbance. The wells were cased with 3 inch 
inside diameter PVC casing to minimize the annular space 
between the casing and the in-hole tools. All five wells 
were cased to 50 ft below land surface, or an elevation of 
about 42 ft. This places the bottom of each well several 
feet below the top of the persistent lower permeability 
layer that forms the bottom of the upper surficial aquifer 



and apparently restricts vertical movement of DNAPL As 
shown in Fig. 9, four wells, identified as 10, 12, 14 and 
16, were completed within the heavily contaminated area 
1 and a "clean" well, 18, was completed at the clean loca- 
tion 5. 

[0052] Each well was logged four times with an active gamma 
(density) tool, a natural gamma tool, and a neutron 
(porosity) logging tool. The active gamma tool utilizes two 
detectors to allow correction for near-hole disturbance. 
The effect of the PVC casing was removed during data 
processing by calibrating the tool with a piece of the cas- 
ing used to construct the wells. The vertical resolution of 
the active gamma tool is about 3-8 inches. 

[0053] The tool for measuring the dielectric constant was 

mounted on a cone-penetrometer tool string (CPT) and 
pushed into the ground with a hydraulic ram. The CPT lo- 
cations were placed as close to the wells as possible. The 
CPT locations are about 1-3 ft from the wells. Each well 
was surrounded by four CPT locations. The CPT locations 
are exemplified by element 20 shown in Fig. 9. The CPT 
tool recorded tip pressure, sleeve resistance, pore pres- 
sure, electrical resistivity, and dielectric constant. 

[0054] The data from the CPT and geophysical logs were interpo- 



lated to consistent 0.1 ft vertical intervals, as variations in 
logging speeds created non-uniform vertical sampling in- 
tervals. A simple linear interpolation was made between 
data points and density and dielectric values were selected 
at corresponding depths at 0.1 ft intervals to calculate 
DNAPL volume. 

[0055] jhe natural gamma logs were used to determine the site 
stratigraphy, as silts and clays are stronger gamma emit- 
ters than clean sands. Using the average of the four log- 
ging runs from well 16 as an example, the gamma stratig- 
raphy suggests a general increase in silt/clay content with 
depth, and three prominent fine-grained layers at eleva- 
tions of about 8 to -10 ft, -18 to 22 ft, and 29 to 33 ft, as 
shown in Fig. 10. The estimated tops of the layers are in- 
dicated on Fig. 10 by dashed lines and the centers of the 
fine-grained layers are indicated on Fig. 10 by arrows 
identifying Unit 1 at 9 ft. Unit 2 at 21 ft, and Unit 3 at 33 
ft. Natural gamma logs from the other four wells show the 
same three principal fine-grained layers. 

[0056] The values measured in the field were the bulk dielectric 
constant, e^, and the bulk density. Porosity was calculated 
from bulk density by assuming 100% water saturation and 
a matrix density of 2.65 g/cc. Plots of e versus elevation 



as shown in Fig. 11 and porosity versus elevation as 
sliown in Fig. 12 for well 16 indicate that both the dielec- 
tric constant and porosity vary significantly with depth. 
Bulk dielectric constant values range from about 14 to 24, 
while porosity varies from about 0.3 to 0.55. A large, low 
value anomaly is apparent in both the porosity and e data 

b 

at about 11 ft. This corresponds to a very low natural 
gamma response at 11 ft, and is correlated with a clean, 
coarse sand. A high porosity, high response occurs at 
about 33 ft. This anomaly corresponds to the highest nat- 
ural gamma counts in well 16, and correlates with a silty, 
clayey fine sand. 
[0057] The neutron log for well 16, as shown in Fig. 13, indicates 
a general decrease in counts per second with depth. There 
are several elevations where the CPS values are lower than 
elevations above or below. These are indicated on Fig. 13 
by the arrows, and occur at 1 ft, -9 ft, -21 ft, -29 ft, and 
33 ft. 

[0058] The presence of DNAPL in pore spaces below the water ta- 
ble reduces the observed bulk dielectric constant from the 
value that would be obtained with no DNAPL present, for 
an observed porosity value. The value of the dielectric 
constant expected for a measured porosity value and 



100% water saturation can be calculated from the equa- 
tion shown in Fig. 6. The vertical displacement of the ob- 
served value of the dielectric constant from the calculated 
value is the result of the presence of DNAPL. Fig. 14 
shows the dielectric constant calculated from Roth's equa- 
tion, assuming 100% water saturation, from field values of 
porosity, and field measured values of the bulk dielectric 
constant. Porosity and dielectric data pairs are from the 
same elevations. 
[0059] The difference between the calculated and observed val- 
ues of the dielectric constant, e andc , is a qualitative 

b t 

measure of the presence of DNAPL. A plot of - versus 
elevation, as shown in Fig. 15, gives a qualitative measure 
of the presence of DNAPL in the subsurface. The plot of 
calculated minus observed dielectric constant suggests 
three principal concentrations of residual DNAPL at eleva- 
tions of 8 ft, -18 ft, and 30 ft. These elevations corre- 
spond to the estimated elevations of the tops of the lower 
permeability layers indicated by the natural gamma log at 
elevations of 9 ft, -21 ft, and 33 ft. These high natural 
gamma response layers have apparently caused the 
DNAPL to accumulate just above them. 
[0060] The observed values of porosity and e at each depth can 

b 



be used to calculate the fractional DNAPL volume, c , us- 

D 

ing the equation of Fig. 6. This calculation gives a quanti- 
tative measure of the volume of DNAPL present as a frac- 
tion of the total volume as shown in Fig. 16. Again, the 
plot indicates three principal concentrations of DNAPL at 
elevations of 8.5 ft, -18 ft, and 30 ft. Fractional DNAPL 
volume at these three elevations is >0.20. DNAPL volumes 
are generally >0.10 between elevations of 5 to 35 ft. 
[0061] The DNAPL saturation is the relative volume of the avail- 
able pore space occupied by DNAPL. It can be easily cal- 
culated by dividing the fractional DNAPL volume, e^, by 
the porosity, 6^. Multiplying the ratio by 100 expresses 
DNAPL saturation as a percentage of available pore space. 
A plot of DNAPL saturation versus depth as shown in Fig. 

17, shows the same three principal DNAPL peaks at 9, - 

18, and 30 ft. At these depths, DNAPL saturation is about 
40-46%. Saturation is generally greater than 25% between 
elevations of 5 to 35 ft. 

[0062] At any given depth, the range in fractional DNAPL is about 
+- 0.05. The most serious source of noise is that the di- 
electric constant and active gamma tools did not pass 
through the same sediment column, as the CPT locations 
were a few feet away from the wells where the gamma 



logs were run. A CPT tool with both an active gamma in- 
strument and the dielectric constant instrument would 
serve to eliminate much of the noise in the data, and 
would make vertical correlation of the dielectric constant 
and porosity measurements more accurate and precise. 
[0063] As illustrated by the flow diagram of Fig. 18, the present 
invention provides a method for the detection and quan- 
tification of a non-aqueous phase liquid that includes 
providing a soil sample 25, measuring the bulk dielectric 
constant of the sample 30, measuring the bulk density of 
the sample 35, estimating the bulk porosity of the sample 
from the measured bulk density 40 and calculating the 
non-aqueous phase liquid saturation from the measured 
bulk dielectric constant and the estimated bulk porosity 
45. 

[0064] It will be seen that the objects set forth above, and those 
made apparent from the foregoing description, are effi- 
ciently attained and since certain changes may be made in 
the above construction without departing from the scope 
of the invention, it is intended that all matters contained 
in the foregoing description or shown in the accompany- 
ing drawings shall be interpreted as illustrative and not in 
a limiting sense. 



[0065] It is also to be understood that the following claims are 
intended to cover all of the generic and specific features 
of the invention herein described, and all statements of 
the scope of the invention which, as a matter of language, 
might be said to fall therebetween. Now that the invention 
has been described. 



